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Filtering,  Smoothing,  and 
Extrapolations  in 
Dose-Response  Experiments: 
With  Application  to  Data  on 
Respiratory  Tumor  in  Rats 

Nozer  D.  Singpurwalla  Jingxian  Chen  ^ 

Abstract 

A  method  for  inference  and  extrapolations  in  certain  dose- 
response.  damage-assessments  and  accelerated  life-testing  stud¬ 
ies  as  been  proposed  by  Meinhold  and  Singpurwalla  in  1986. 

The  method  is  based  on  a  use  of  the  Kalman-filter  algorithm 
and  involves  the  double  lognormal  as  the  distributional  assump¬ 
tion.  In  this  paper,  we  discuss  issues  pertaining  to  a  practical 
implementation  of  this  methodology.  This  involves  some  in¬ 
sights  based  on  a  simulation  study  about  the  specification  of 
prior  parameters  and  an  application  of  the  proposed  methodol¬ 
ogy  to  some  published  data  on  doses  of  bischloromethyl  ether 
administered  to  rats. 

Key  Words  and  Phrases:  Accelerated  Testing,  Bioassay, 

Kalman  Filtering,  Bayesian  Methodology,  Damage  Assessment. 

1,  INTRODUCTION  AND  OVERVIEW 

Let  X  denote  a  dose  or  a  stress  that  is  applied  to  a  biological  or  an  engi¬ 
neering  system,  and  suppose  that  x  takes  values  in  [0,oo).  Let  ^(x)  be  the 
response  to  x.  and  suppose  that  0  <  y^(x)  <  1;  V  (x)  could  be  viewed  as  the 
propensity  of  each  item  to  respond  to  x,  or  the  extent  of  damage  incurred 
by  each  item  when  subjected  to  stress  x.  Thus,  for  example,  Y{x)  =  1(0) 
could  imply  total  resistance  (demolition)  to  (under)  x.  Often,  it  is  true  that 
we  are  able  to  test  more  than  one  item  at  any  x,  but  that  it  may  not  be 
possible  to  repeat  the  test  at  any  x,  because  doses  and  stresses  cannot  be 
exactly  controlled.  In  what  follows,  we  shall  assume  that  the  values  of  x 
are  known  to  an  experimenter,  but  that  the  Y{x)  are  unknown  “states  of 
nature”  about  which  it  is  desired  to  make  inferences.  Such  inferences  are 
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based  on  background  knowledge  about  a  particular  scenario  under  consid¬ 
eration,  an  assumed  relationship  between  and  x,  and  tests  conducted 

at  several  values  of  x.  The  tests  conducted  at  the  several  values  of  x  give 
us  information  about  Y{x),  but  that  Y[x)  cannot  be  directly  observed. 

It  is  often  the  case  that  Y{x)  responds  to  2:  in  a  nonlinear  fashion,  and 
in  view  of  this,  plus  considerations  of  the  type  cited  in  Section  1.1,  we 
propose,  as  a  model  for  relating  Y{x)  and  x,  the  relationship 

Y{x)  =  exp[— cv(i),  ^(2;)  >  0;  (1-1) 

«(x)  and  3(x)  are  unknown  parameters  which  depend  on  x.  In  what  follows, 
£’(Y(x))  will  denote  expectation  of  the  unknown  quantity  V'(x). 

1.1  Arguments  Supporting  Choice  of  the  Relationship 
As  a  special  case  of  (1.1),  suppose  that 

a(x)  =  rt  and  0(x)  =  6,  for  all  x  >  0. 

Then 

y'(x)  =  exp(— ftx'^),  a,  >  0;  (1-2) 

the  right  hand  side  of  the  above  is  the  survival  function  of  a  Weibull  distri¬ 
bution.  A  virtue  of  (1.2)  is  flexibility  in  expressing  a  wide  class  of  subjective 
opinions  about  the  dose-response  and  stress-damage  relationships.  For  ex¬ 
ample,  it  has  been  recommended  for  use  in  food  safety  assessment  studies 
and  clinical  trials  in  the  biological  scenarios  [cf.  Final  Report  of  the  SCF- 
SC  (1980)],  and  in  studying  the  effects  of  underwater  nuclear  explosions 
on  submarine  miniatures  in  the  engineering  scenario  [cf.  McDonald  (1989), 
Shaked  and  Singpurwalla  (1990)].  The  relationship  (1.2)  [and  also  (1.1)] 
can  be  linearized  enabling  one  to  employ  standard  filtering  techniques,  and 
when  used  in  connection  with  data  sets  on  human  risk  assessment  studies, 
it  gives  an  estimated  risk  at  low  doses  which  lies  between  the  estimates  for 
the  “gamma  multi-hit”  and  the  “.Armitage-Doll”  models.  Also,  the  rela¬ 
tionship  (1.2)  has  been  reported  to  give  good  empirical  fits  to  the  data.  In 
any  particular  application,  it  is  suggested  that  the  statistical  analyst  and 
the  subject  matter  specialist  examine  plots  of  the  Weibull  dose-response 
curves  for  several  combinations  of  values  of  a  and  0,  and  choose  that  curve 
or  curves  that  best  describe  their  judgments  of  the  dose-response  relation¬ 
ship.  The  above  can  be  most  effectively  done  on  personal  computers. 

1.2  Statement  of  the  Problem 

Suppose  that  at  T  distinct  dose  levels  xj  >  X2  >  ■  •  ■  >  2;t  an  item  or  sev¬ 
eral  items  are  tested  and  the  corresponding  observed  responses  y{xi),  y(x2), 
. .  .,y{xT),  recorded;  note  that  t/(x,)  provides  information  about  the  unob¬ 
servable  Y{xi),i  =  I, . .  .,T.  We  are  required  to; 


Dose  Response  Experiments 


279 


i)  Make  statements  of  uncertainty  about  the  true  responses  V'(a:),  at  any 

X,  including  the  i.’s,  i=  at  which  the  tests  are  conducted. 

ii)  Assess  y'{a;T+i)>  the  response  at  a  very  low  dose  ir+i.  where  ir+i 
<<  XT 

1.3  The  Proposed  Approach 

Our  approach  for  addressing  the  above  issues  is  Bayesian,  and  involves  a 
use  of  the  technology  of  Kalman  filtering  and  Kalman  filter  smoothing.  The 
underpinnings  of  our  approach  are  outlined  in  Meinhold  and  Singpurwalla 
(1987),  where  due  to  limitations  of  space  and  the  need  for  an  exposito¬ 
ry  focus,  applications  of  the  methodology  to  practical  scenarios  was  not 
undertaken.  An  aim  of  this  paper  is  to  fill  the  above  gap  and  to  attempt 
to  demonstrate  the  potential  usefulness  of  the  Kalman  filter  approach  to 
problems  of  the  kind  considered  here.  It  is  helpful  to  point  out  that  Black- 
well  and  Singpurwalla  (1988)  undertake  an  endeavor  similar  to  the  one 
discussed  here,  but  focus  attention  on  the  case  of  exponentially  distributed 
lifetimes:  furthermore,  they  do  not  address  the  issue  of  smoothing  that  is 
relevant  to  problems  of  this  type. 

1.4  Overview  of  Paper 

In  Section  2  we  motivate  the  Kalman  filter  model,  state  the  distribution¬ 
al  assumptions,  and  present  the  necessary  smoothing  and  extrapolation 
formulae.  In  Section  .3  we  describe  approaches  for  obtaining  the  starting 
values  of  the  Kalman  filter  algorithm  (i.e.  specifying  the  prior  parameters), 
and  in  Section  4  we  describe  an  application  of  our  approach  to  some  data 
on  respiratory  tumor  of  rats  subjected  to  doses  of  bischloromethyl  ether.  In 
Section  5  we  offer  some  conclusions  and  suggest  some  direction  for  future 
research  along  the  lines  outlined  here. 

2.  THE  FILTERING  AND  SMOOTHING  MODEL 

Our  review  of  the  literature  suggests  that  there  is  a  dearth  of  dose- 
response  relationships  based  on  pharmacokinetic,  oncological  or  engineer¬ 
ing  considerations.  Thus  one’s  choice  of  (1.2),  the  Weibull  survival  function, 
is  at  best  a  reasonable  approximation,  and  hence  one  should  incorporate 
into  one’s  analyses  some  meeisure  of  uncertainty  about  such  approxima¬ 
tions.  Furthermore,  there  is  no  recison  to  sacrifice  flexibility  by  making  the 
Weibull  scale  parameter  cv,  and  the  shape  parameter  /?,  invariant  with  the 
dose  X.  On  the  contrary,  there  is  evidence  in  the  engineering  sciences,  that 
high  stresses  cause  a  change  in  the  basic  failure  mechanisms,  making  it 
reasonable  to  cissume  that  a  and  /?  depend  on  x  -  thus  our  choice  of  the 
relationship  (1.1).  The  dynamic  nature  of  a  and  0  introduces  the  novelty  of 
our  approach  over  the  currently  used  ones.  Once  the  above  is  undertaken, 
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our  set-up  fits  into  the  general  form  of  a  Kalman-filter  model  -  provided 
that  some  details  which  facilitate  an  iterative  computational  scheme  are 
attended  to. 


2.1  The  Kalman- Filter  Model  and  Distributional  Considerations 

Assume  that  E{y{x))  =  exp(-o;(a:)x‘^*^)),  and  suppose  that  we  require 
d^f 

y'{x)  =  log{— log y(z)},  to  be  such  that 

y’{x)  ~  .\f{n{x),a~{x)), 

where  denotes  “normally  distributed  with  mean  (variance) 

fi{x){cr~{x))” .  Then,  it  can  be  seen  that  y{x)  must  have  a  “double  lognor¬ 
mal  distribution”  [cf.  Meinhold  and  Singpurwalla  (1987)]  with  parameters 
/i(x)  and  o’^(x),  where  /i(x)  is  the  median  of  the  distribution  of  y{x).  The 
distribution  cited  above,  has  a  density  which  is  flexible  enough  to  express 
a  wide  variety  of  subjective  opinions  about  y(x),  and  is  an  attractive  alter¬ 
native  to  the  beta  density  for  modelling  data  on  proportions  -  see  Ahsan- 
ullah  and  Holland  (1989).  Motivated  by  the  fact  that  when  (T“(x)  is  small. 
E{y{x))  Si  exp{— exp(/i(x))},  we  propose,  as  the  observation  equation  of 
the  Kalman  filter 


y'{x)  =  (l.logx) 


+  i^(i), 


(2.1) 


where  i/(x)  ~  .'V'(0,  V^(x)),  V(x)  is  the  conditional  variance  of  y*  (x),.  (7, /?)^ 
=  (7(a:)>/?(ic))',  and  7(x)  =  loga(x). 

The  set-up  (2.1)  also  implies  that  y(x)  =  (V'(x))''(^\  where  logA(x)  = 
i^(x),  and  2/’(x)  =  log{— log 2/(x)};  the  innovation  A(x)  has  a  lognormal 
distribution  with  parameters  0  and  cr-{x).  For  the  system  equation  of  the 
Kalman  filter,  we  propose  the  “steady  model”  [cf.  Meinhold  and  Singpur¬ 
walla  (1983)] 


’7' 

J. 

(r) 

J. 

+  w{x), 


(2.2) 


where  to(x)  ~  A/”(0,  M^(x)),  and  W(x)  is  the  variance-covariance  matrix 
of  7(x)  and  /?(x)  conditional  on  7(x  —  1)  and  0(x  —  1);  //(x)  is  assumed 
independent  of  tn(x);  (x—  1)  is  the  dose  previous  to  dose  x.  The  values  V{x) 
and  the  four  entries  of  W(x)  must  be  specified  by  a  user.  Also  to  be  specified 
by  the  user  are  the  “starting  values”  of  the  Kalman  filter  algorithm;  these 
are  7(0),  /?(0),  V{xo)  and  W(xo).  A  strategy  for  specifying  these  is  outlined 
in  Section  3. 
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2.2  Results  from  Filtering  and  Smoothing 

For  compactness  of  notation,  we  let  y^{xj)  =  2/| ,  [7./^]ry  =  = 

Fj,  V{xj)  =  Vj,  W{xj)  =  Wj,  ty(xj)  =  Uj  and  w(xj)  =  iVj.  Then  (2.1) 
and  (2.2)  can  be  written  as 

y]  =  FjOj  FUj,  and 
Oj  —  Gj6j-  \  +Wj, 

where  Gj  is  a  2  x  2  identity  matrix.  If  we  assume  that  dj-i  ~  jV{9j.u 
Sj-i),  then,  upon  observing  y‘ ,  dj  ~  M  (dj,  Ej),  and 

6j  =  Gjdj-i  +  Kj{y*  —  FjGj6j-i),  Sj  =  ( J  —  KjFj)Rj,  (2.3) 
where 

Rj  —  Gj  Ej-i  G'j  +  W  j,  Kj  =  Rj  Fj  {Fj  Rj  Fj  +  V’j )  , 

and  I  is  an  identity  matrix.  The  relationships  in  (2.3)  are  referred  to  as  the 
forward-recuTrence  equations  of  the  Kalman  filter.  Note  that,  in  the  above 
scheme,  inference  for  6j  is  based  upon  the  previous  data,  yj ,  yj^\,  ■  ■  ■  ,yi, 
j  =  1,2,,,.,T,  only.  Thus,  with  the  exception  of  Or,  inference  for  the 
other  Oj’s  is  not  based  on  all  the  data.  Should  we  desire  inference  for 
Bj  based  on  aU  the  data  then  we  need  to  smooth  using  the 

backward-recurrence  equations  [see  Appendix  A  of  Meinhold  and  Singpur- 
walla  (1987)],  whereby 

r  OjiT)  =  djij)  +  Jj+,[9j+iiT)  -  Gj+i9j{j)],  (2  4) 

\  Ej{T)  =  Sj{j)  -  Jj+i[Ej+i{T)  -  Ri+i]J'j+w 

where  Jj  =  Ej-i{j-\)  G'j  RJ^  and  0,(l;)[-5^j(^')]  is  the  mean  [covariance] 
of  the  normal  distribution  of  9j  based  on  y[ , . . . ,  y’J.  The  standard  Kalman- 
filter  algorithm  enables  us  to  produce  ( 1),  ^2(2),  •  •  •  i  9t{T)  in  an  efficient 
manner;  see  (2.3). 

Once  9t{T)  is  obtained,  inference  about  T(it+i)  follows  from  the  fact 
that  log(-log(y'(2:T+i)))  ~2V’(Ft+i  9t{T),  Fr+iEriT)  Fj+J,  and  so 

F(y(a;T+i))  ~  exp(-exp(FT+i^T(T))).  (2.5) 

3.  SPECIFICATION  OF  INPUTS  FOR  THE  KALMAN  FILTER 

Let  xo  and  x_i,  x_i  >  xq  >  xi,  be  two  dose  levels  at  which  the  subject 
matter  specialist  heis  the  most  knowledge  about  the  responses,  y(xo)  and 
y(x_i)  respectively.  Typically,  these  would  be  very  large  dose  levels;  also, 
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we  require  that  xq  and  a:_i  be  close  to  each  other.  Suppose  that,  in  the 
opinion  of  the  specialist,  the  most  likely  values  of  ^(so)  and  y'(x_i)  are 
V'(2o)  =1—^0  and  Y (i_  i )  =  1  — i5i ,  where  6o  and  6i  are  large  with  i5o  <  6\ ; 
what  we  have  in  mind  is  6o  =  .995  and  (5i  =  .999.  Then,  for  xq  close  to 
I-  1 ,  we  may  write 


log{-log(l  -5o)}  ss7(a:o)  +  /?(a:o)logio, 


log{-log(l  -^i)}  «7(a;o)  +  /?(a;o)loga:_i, 
from  which  it  follows  that 


(3.1) 


.  _  {10g(-log(l  -(Sq))  -  log(-log(l  -  ^;))} 

°  {logxo  -  logx_i} 

and 

'/(xo)  =  log(-log(l  -5o))  -/f(xo)logxo. 

The  values  l3(xo)  and  ‘/(xq)  thus  computed  will  become  the  starting  values 
for  the  Kalman-filter  algorithm. 

Our  next  task  is  to  pin  down  V^(xo), . . . ,  V(x7’),  the  variances  of  the 
innovations  r'fx,),  i  =  0.  .  .  .  ,  T.  In  specifying  the  above,  two  considerations 
must  be  borne  in  mind.  The  first,  is  that  since  y{x)  is  between  0  and  1, 
the  variance  of  A(x)  must  be  the  smallest  when  V'(x)  is  either  0  or  1,  and 
the  largest  when  y(x)  =  .5.  Second,  under  binomial  testing,  the  variance 
of  y(x)  is  approximately  of  the  form  m(x)  (1  —  m(x))/n(x),  where  m(x) 
is  the  median  of  the  distribution  of  y(x)  and  n(x)  is  the  number  of  units 
tested  at  dose  x.  Following  a  line  of  reasoning  given  in  Appendix  A,  a  first 
order  approximation  to  V(x)  =  <x'(x)  is  of  the  form 

<7-{x)  K  log  1^(1  +  [1  +4(1  -  m(x))/7i(x)m(x)(log(m(x)))-]2)|  .  (3.2) 


Note  that  in  the  above  expression,  cr-(x)  —  0  as  m{x)  —  0[1]  provided 
that  ri(x)  >  (m(x))"'  [(1  —  ni(x))~'],  and  so  as  is  commonly  done  by  en¬ 
gineers,  we  may  set  V(xo)  =  0.  In  (3.2)  above,  m(x)  =  exp{- e.xp(/i(x))}, 
with  /i(x)  replaced  by  its  predicted  value  /i(x)  obtained  from  the  Kalman- 
filter  algorithm.  Specifically, 


+(2^i)  =  7(^0)  + /?(xo)  log xi,  and 

fi{xi)  =  FiGi,  i=l,...,T, 

where  di  is  given  by  (2.3),  the  forward  mechanism  of  the  Kalman  filter. 

Our  next  step  is  to  pin  down  W(xo),  the  variance  covariance  matrix  of 
the  starting  values  7(xo)  and  P(xo).  For  this,  we  require  that  the  subject 
matter  specialist  specify,  in  addition  to  Y (xq)  and  Y (x_i),  i/i,  and  p,  the 
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variances  of  y(a:o),  y(i-i)  and  the  correlation  between  y^;^,  and  y 
as  perceived  by  the  specialist.  The  variances  and  1^2  will  typically  be 
very  small,  something  like  (.001)"  and  (.0003)",  respectively,  and  p  should 
be  nonnegative,  say  something  of  the  order  of  .7  or  .8.  Once  the  above  are 
done,  the  relationship  (3.1)  can  be  used  to  show  that  if  f>  =  i  and 

if  Z  =  [Var(7(ro)),  Var(/d(xo)),  Cov(7(xo), /3{xo))]',  then  Z  =  A-‘b,  where 
the  matrix  A  is  of  the  form 


i'i  +  i'2  +  2pv 


(logxo)^ 

(logx_i)- 

pflogio+logr-i) 
t' 1  +  2  +  2  p 


2  log  XO 
21ogx_i 

4p(lopCiro+loK3:^i) 
i  +  *'2  +  2 P \/ 


Once  Z  is  known,  the  matrix  I7(xo)  can  be  constructed,  and  this  en¬ 
ables  us  to  undertake  the  first  iteration  of  the  Kalman-filter  algorithm. 
Subsequent  iterations  of  the  algorithm  require  that  we  specify  W(z;), 
j  =  1. . . . ,  T.  For  this,  we  propose,  based  on  some  simulation  studies,  that 

A 

W(x^)  =  ..5(1.5)^(xj-i  -  Xj)Sixo).  (3.3) 


The  multipliers  of  S(xo)  given  above,  reflect  the  following  considerations; 

i)  The  factor  .5  (or  any  other  number  less  than  1)  reflects  the  fact  that 
our  uncertainty  about  the  parameters,  subsequent  to  observing  data, 
should  be  less  than  our  uncertainty  about  them  prior  to  the  data. 


ii)  The  factor  (1.5)-'  (or  for  that  matter  any  other  number  greater  than 
1)  reflects  the  fact  that  our  uncertainty  about  the  parameters  should 
increase  as  we  get  closer  to  the  low  dose  levels. 


iii)  The  factor  {xj-i-Xj)  reflects  the  fact  that  if  the  separation  between 
two  consecutive  doses  is  large  then  our  uncertainty  about  the  model 
parameters  should  increaise. 

It  should  be  clear  from  the  above,  that  the  scheme  proposed  here  could  be 
automated  once  y(a:o),  1^1,  y(x_i),  1^2  and  p  are  specified. 


4.  APPLICATION-Data  011  Respiratory  Tumor  in  Rats 


For  purposes  of  illustration,  we  consider  some  data  on  doses  (the  num¬ 
ber  of  six  hour  exposures  by  inhalation  of  100  parts  per  billion)  of  bis- 
chloromethyl  ether  administered  to  rats.  This  data  has  been  abstracted 
from  the  Final  Report  of  the  SCFSC  (1980),  wherein  its  analysis  using  the 
Weibiill  survival  model  htis  been  advocated  and  undertaken.  In  Table  4.1 
below,  we  present  our  data;  the  response  is  0  if  a  rat  develops  respiratory 
tumor  and  1  otherwise  -  thus  43/46  denotes  the  fact  that  43  out  of  46 
rats  have  not  developed  tumors.  Our  aim  is  to  predict  the  responses  at  low 
doses,  say  7,  5,  2  and  1. 
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Table  4.1.  Dose-Response  Data  on  Respiratory  Tumors  in  Rats 


Dose:  x,- 

10 

20 

40 

60 

80 

100 

Observed 

Response 

40/41 

j 

43/46 

14/18 

14/18 

19/34 

8/20 

We  let  xi  =  100,^2  =  80 . xe  =  10,  and  y(xi)  =  40/41, 1/(12)  = 

43/46, . . . ,  y(xe)  =  8/20.  For  xq  and  i_i  we  choose  the  doses  250  and  280 
respectively,  and  for  ^0  and  (5i  we  choose  .995  and  .999,  respectively.  Our 
other  choices  follow  the  recommendations  of  Section  3;  specifically,  1/1  = 
(.001)“,  1^2  =  (.0003)“  and  p  =  .7.  Substituting  the  above  in  the  formulae 
for  /i(xo),  '/(xo)  and  A,  we  obtain,  as  starting  values  for  the  Kalman  filter, 
the  following; 

7(xo)  =  -11.256,/5(xo)  =  2.3405,  Var(7(xo))  =  1.668  x  10-^, 

Var(/3(xo))  =  5.217  x  10~*  and  Cov  {j'(xQ),  /)(xq))  =  —2.95  x  10“"*. 

In  Table  4.2  below,  we  give  filtered  and  smoothed  estimates  of  the  pa¬ 
rameters  7(x,)  and  d(x,),  2  =  1, . .  .,6.  The  filtered  estimates  are  obtained 
via  the  forward  recurrence  equations  (2.3)  whereas  the  smoothed  estimates 
are  obtained  via  the  backward  recurrence  equations  (2.4).  The  entries  in 
Table  4.2  indicate  that  smoothing  does  have  an  effect  on  the  filtered  esti¬ 
mates  and  that  the  estimates  of  7(x,  )  and  0(xi)  do  change  with  x,-.  Note 
that  when  i  =  6.  that  is,  for  Xis  =  10,  the  smoothed  estimate  is  indeed  the 
filtered  estimate. 


Table  4.2.  Filtered  and  Smoothed  Estimates  of  7(xi)  and  /?(x,) 


Dose:  x,- 

100 

80 

60 

40 

20 

10 

Filtered  Est. 
of  7(x,) 

-10.652 

-9.563 

-10.079 

-7.702 

-7.965 

-7.604 

Smoothed  Est. 
of  7(xi) 

-9.801 

-9.780 

-9.651 

-7.704 

-7.965 

-7.604 

Filtered  Est. 
of  3{xi) 

2.2.34 

2.041 

2.132 

1.712 

1.759 

1.695 

Smoothed  Est. 
of  3{xi) 

2.083 

2.080 

2.057 

1.713 

1.759 

1.695 

In  Table  4.3  we  give  the  observed  values  j/(x,)  and  the  predicted,  filtered 
and  smoothed  values  of  T (x,),  z  =  1 , . . . ,  6.  Note  that 

i)  £'(y(x,)|j/(x,_i ),...,  y(xi )),  is  the  predicted  value  of  y(x,),  where 
y(xo)  =  y(xo), 

ii)  E(V (xi)jy(xi), . . . .,  y(xo))  is  the  filtered  value  of  Y (x,),  and 
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iii)  E{Y[xi)\y{x^), . yi^o))  is  the  smoothed  value  of  y(ii)- 

The  above  conditional  expectations  can  be  obtained  via  formulae  analogous 
to  (2.5). 

Table  4.3.  Predicted,  Filtered,  Smoothed  and  Observed  Values  of  Y{xi) 


Dose:  x,- 

100 

SO 

60 

40 

20 

10 

Observed  Val. 

.4 

.5588 

.7778 

.7778 

.9348 

.9756 

of  y{xi) 
Predicted  Val. 

.5377 

.6559 

.7411 

.8963 

.9265 

.9756 

of  y(x.  ) 

Filtered  Val. 

.4995 

.5834 

.7712 

.7788 

.9348 

.9756 

ofy(x.) 
Smoothed  Val. 
of  y  (x,  ) 

.4439 

.5986 

.7465 

.7789 

.9348 

.9756 

We  note  from  Table  4.3,  that  smoothing  has  a  tendency  to  lower  the 
filtered  values  of  Y (xi)  and  that  the  predicted  values  have  a  tendency  to 
be  larger  than  the  observed  values.  Finally,  in  Table  4.4,  we  give  predicted 
values  of  Y{xi),  for  x,-  =  7,  5,  2  and  1.  low  doses,  at  which  no  testing  was 
done  and  at  which  inference  is  most  crucial.  Note  that  here  the  predicted 
values  y{xi)  are  given  by  E(y (x’,)|!/(x6)-  ■  •  •  ■  '/(■co)))  and  that  this  quantity 
is  obtained  via  (2.5).  Also  given  in  Table  4.4  are  the  90%  Probability  of 
Coverage  Intervals  (L,',L't)  for  the  above  V''(x,)’s. 


Table  4.4.  Predicted  Values  of  y{x,)  at  Low  Doses 


Dose:  x,- 

7 

5 

2 

1 

y{xi) 

.9866 

.9924 

90%  PCI 

L 

.98658 

.9984 

.9995 

for  y(x,) 

U 

.98662 

.9984 

.9995 

4.1.  Comparison  with  Maximum  Likelihood  Approach 

It  is  of  interest  to  compare  the  results  produced  by  our  approach  with 
those  produced  via  the  conventional  approach,  in  which  the  relationship 
(1.2)  is  assumed,  and  a  and  0  estimated  by  the  method  of  maximum  like¬ 
lihood.  It  can  be  easily  verified  that  a  and  0,  the  maximum  likelihood 
estimates  of  cv  and  0,  respectively,  are  d  =  7.496  x  10~'*,  and  0  =  1.513. 
Replacing  n  and  0  hy  a  and  0  respectively,  in  (11),  we  obtain  the  max¬ 
imum  likelihood  estimates  of  y(x,).  These  are  shown  in  Table  4.5;  also 
shown  there,  for  purposes  of  comparison,  are  the  smoothed  values  of  y(x,) 
and  the  observed  values  y(x,).  We  note  from  Table  4.5,  that  whereas  the 
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differences  between  the  maximum  likelihood  estimates  of  Y{xi)  and  the  s- 
moothed  values  of  y(x,')  are  significant  at  the  higher  doses,  the  differences 
at  the  lower  doses,  particularly  the  doses  at  which  predictions  are  sought, 
are  negligible  if  not  nonexistent.  Thus  it  appears  that  a  use  of  our  approach 
does  not  lead  one  to  conclusions  that  would  be  significantly  different  from 
those  obtained  via  a  more  conventional  approach,  except  that  should  the 
situation  so  demand,  our  approach  would  provide  a  greater  flexibility  than 
the  conventional  one.  Finally,  we  also  note  that  the  smoothed  values  of 
V'(x,  )  tend  to  be  larger  than  those  of  the  maximum  likelihood  estimates. 


Table  4. .5.  A  Comparison  of  the  Observed  and  Smoothed  Values  and 
Maximum  Likelihood  Estimates  of  Y(xi) 


Dose:  x,- 

Observed  Value 
of  y(xi) 

.Smoothed  Value 
of  yfx,) 

Max.  Likelihood  Est. 
o{Y(xi) 

100 

.4000 

.4439 

.4515 

80 

,.5.588 

.5986 

.5667 

60 

.7778 

.7465 

.6926 

40 

.7778 

.7789 

.8197 

20 

.9348 

.9348 

.9327 

10 

.9756 

.9756 

.9759 

7 

.9866 

.9859 

5 

.9924 

.9915 

2 

.9984 

.9979 

1 

.9995 

.9993 

5.  SUMMARY  AND  CONCLUSIONS 

It  appears  to  us  that  using  the  dynamic  linear  model  set-up  of  Kalman 
filtering  is  a  potentially  useful  approach  for  making  inference  under  dose- 
response  experiments.  Its  chief  virtues  are  flexibility  in  modelling  -  specifi¬ 
cally,  making  the  parameters  dose  dependent  -  and  the  ease  of  undertaking 
inference.  Its  chief  disadvantage,  especially  to  a  non-Bayesian.  is  the  need  to 
specify  the  starting  values.  However,  as  indicated  in  Section  3,  some  general 
guidelines  can  be  followed,  and  once  this  is  done  the  procedure  is  almost 
automatic.  A  computer  code  which  facilitates  the  required  computations 
are  given  by  Chen  and  Campodonico  (1989).  An  issue  that  remains  to  be 
addressed,  and  one  that  we  have  not  been  able  to  satisfactorily  undertake, 
is  that  pertaining  to  the  determination  of  a  “safe  dose"’.  Specifically,  what 
is  needed  is  inference  about  x,  when  Y{x,)  is  specified;  Y{x,)  is  typically, 
a  number  close  to  1,  say  .999  or  .9999.  Finally,  regarding  further  research 
along  the  above  lines,  it  would  be  desirable  to  develop  an  inference  mech¬ 
anism  which  does  not  rely,  as  heavily  as  we  have,  on  a  use  of  the  Gaussian 
distribution. 
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APPENDIX  A 

THE  SPECIFICATION  OF  a-(x) 

Let  Z{x)  =  -logy(z)  =  (- log  F(a:))A(a:),  then  Z(x)  has  a  lognormal 
distribution  with  parameters  log(— log  F(x))  and  (t-{x). 

Since  y(x)  =  m(x)  •  exp{-(Z{x)  +  logm(x))},  we  have  by  a  Taylors 
series  expansion 


y(x)  =  m(x)^(-l)*(Z(x)  +  logm(x))V^! 

Jc=0 

w  m(x)(l  -  (Z(x)  +  logm(x))) , 
from  which  it  follows  that 


Var  (y(x))  «  m-(x)  •  Var  (Z(x)) 

=  r7i^(x)(logy(x))-e‘''<"He"'^"^- !)• 

However,  /(x)  =  mix),  since 

.5  =  Pivix)  <  m(x)) 

=  P(y(x)'^''"'  <  m(x)) 

=  P(A{x)  >  logm(x)/logy(x)). 

Therefore,  log  m(x)/  log  y(x)  equals  to  the  median  of  A(x),  which  is  equal 
to  1.  Now  (3.2)  follows  from  simple  algebra  and  the  fact  that  Var(?/(x))  = 
m(x)(l  -  m(x))/7i(x). 
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